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Abstract 

We present a new representation of spin operators in terms of bosonic 
creation-annihilation operators. This representation allows us to formulate 
a new field-theoretical description of spin systems which is free of any con- 
straints. The corresponding functional integral representations for thermo- 
dynamic quantities are given and the application to investigations of Long 
Range Order in the system is discussed. 
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1 Introduction 

There are two sources of motivation to search for bosonic representations of spin 
systems and systems of truncated oscillators: the first is technical while the sec- 
ond is of principle. Indeed, a general problem of any perturbative investigation of 
spin systems is the complicated diagram technique which originate from spin-spin 
commutation relations. On the other hand, for systems with the Hamiltonian for- 
mulated in terms of the bosonic or fermionic creation -annihilation operators, the 
diagram technique is standard and straightforward. That is why we need a bosonic 
representation for the spin operators to cast the complicated technique into the com- 
mon form and use the field-theoretical machinery. Another set of problems where 
the bosonic treatment is vital is when looking for Long Range Order (LRO) in the 
systems. It is well-known that LRO is reflected in apearance of anomal averages. It 
is always very tempting to reformulate the problem in such a way that the anoma- 
lous averages become amplitudes of a Bose-condensate of some auxiliary bosons. 
This was a guidline, for example, in Ref. [ffl where constraint-free representations 
was found for Paulions to predict the Bose-condensation of Frenkel excitons. In this 
paper we go along a similar line and develop a constraint-free description for arbi- 
trary spin system. To this end we make use the approach developed for truncated 
oscillators in 0. 
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We have to note that there are several transformations that expresse the spin oper- 
ators in terms of the bosonic or fermionic ones |3], f|]. However, all of them require 
either the restriction of the bosonic Hilbert space which leads to constraints for the 
bosonic system or restrict the study to ID systems. The constraints do not cause 
any problem unless the systems are treated exactly. Since it is very difficult to get 
exact results for spin systems, some type of approximations should be used. The 
most popular approximation scheme is based on the mean field description. At this 
point drawbacks of the constrained description emerge. Indeed, the mean field ap- 
proximation does not treat local (on-site) constraints in a proper way. It means that 
instead of many local constraints only one global constraint appears. All together it 
leads to the problems of the account of unphysical local fluctuations. This effectively 
returns us to the local constraints and explains the importance of the constraint-free 
formulation of the mapping from spin systems to bosonic ones. 
Similar to the approach of the sigma-model with Wess-Zumino term ||] we treat 
the constraint on the number of particles on each site exactly. To do this we use the 
mapping of the orthogonal sum of identical copies of the lattice spin space of states 
to the bosonic space of states. In this mapping spin operators are represented in 
the form of a power series of the bosonic creation and annihilation operators. This 
compels us to deal with infinite series of different vertices in the diagram technique. 
The choice of relevant contributions in such series should be dictated as usually by 
features of the concrete problem. 

2 Mapping of spins to bosons without constraint 

In this section we will describe the mapping from the system of lattice spins to the 
auxiliary bosonic system. The goal is to escape the introduction of a constraint. To 
do this we will embed an infinite number of copies of the finite dimensional space 
of states in the bosonic space of states and then proceed with the consideration of 
this new (auxiliary) bosonic space. 

To explain this in detail, let us first of all consider one degree of freedom (i.e. a 
single site). Spin operators obey the following commutation relations (for spin m/2): 

{s-) m+1 







S + S- = 2S Z , {S + ) + = S- , (S + ) m+1 

S + = s x + iS v , S- = S x -tS y . (1) 

Operators S + , S~ and S z have the following matrix form in the m + 1-dimensional 
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with basis {|0), |1), |m)} and the obvious notations. 

Now we introduce the infinite orthogonal sum Ti b = © 7~Ls,n of such m- dimen- 
sional Hilbert spaces Hs, n with basis {{|0), |1), . . . , \m — 1)}, . . . , {\nm + 1), \nm + 
2), . . . , \nm + m)}, . . .}. The extensions of operators S + , S~ and S z in this space 
have the form: 

S + = diag(S + , S + , ...) , S~ = diag{S~ , S~ , ...) , 

S z = diag(S z ,S z ,...) . 

It follows that all thermodynamic quantities calculated with operators S + , S~, S z 
are exactly the same as those calculated with the original operators S + , S~ , S z . 
Indeed, for example, 

g+g ^ Sp(S+S-e-K E -ri s+s ~-) 

coincides with the same expressions but without hats due to the block structure 
of our operators (we should add that the partition functions differ by an infinite 
numerical constant which does not affect observable physical quantities). The con- 
clusion is still valid if we start with a lattice of spins and then introduce hats for 
the operators. 



3 Bosonic representation for spin operators 

Let us now derive the relations for matrix elements of operators S + , S~ and S z . To 
do this we will follow the method proposed by Chernyak in Ref. for Paulions. 
The main point of the method is to use the projection operator on the vacuum state 
of the auxiliary boson system, i.e on the vector |0). This projection operator V has 
the following expression in terms of the bosonic creation and annihilation operators: 

00 f— 1V 

V = J2 ^-j^ib+yb 1 =. exp(-b + b) : . 
i=o lj - 

We now can use this representation to construct the operators S + , S~ and S z which 
obey algebra (|l]). Indeed, it is easy to check from the matrix form that the following 
relations hold: 

co m— 1 \ (k ~\- 1) (tTI k) 

g+ _ Qj+^mn+k+l'pjjTrm+k V ^ '_2 = + 

^ ^ {mn + k)\y/mn + k + 1 ' 



n=0 k=0 
oo m 



n=ofe=o (mn + k)\ 

It is obvious that these relations satisfy the algebra (|l|). For the particular case 
m = 2 our formulae reduce to the formulae originally obtained by Chernyak M for 
the case of paulionic operators. 
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4 Thermodynamics in functional integral repre- 
sentation 

The formulae considered above can be applied to construct the Hamiltonian of the 
auxiliary bosonic system. Let us start with the following Hamiltonian H$ of spins 
on a lattice: 

Hs = X x ikSiSt + X YikSfSk + X z ikStS z k + 

i^fc j^fc j^fc 

Using operators S + , S~ one can cast it in the following form: 

H s = \T,( x *h + yik)StSu+\^{{X ih -Y ik )StSt + h.c.} + 

iy^k ij^k 

+ X Z ikS-S z k + ~ £ {(h ix - ih iy )Sf + h.c.} + x h iz S* . 

iy±-k i i 

The corresponding Hamiltonian of the auxiliary bosons based on the relations (|2]) 
has the form: 

H = lT,(Xk + yk)btS ik b k + ^Y,{( X ^~yk)btbiS tk + h.c.} + 

iy^k ij^k 
oo 

+ X^ X a{l)a{m){bf)\btrb\b^ + 

i^k l,m=0 

1 ( oo ^ oo 

+ 5 x - «^) x + fc.c. + x x «co( w 

/ j l Z=0 ) i 1=0 

Here the following notations have been introduced: 

oo 

S lk = X A(l)A(m)(biY(birb[b^ , 

l,m=0 

A(Z) = v/2 X X V 



fc=0 



I i-fcl 



^ (/ — mn — k)\ (mn + k)\y/mn + k + 1 



(_l)i-mn-fe (_ m/2 + l) 

Using the standard procedure, we can write down the functional integral represen- 
tation of the partition function and correlators of the auxiliary bosonic system and 
the original system of truncated oscillators. For example, according to the definition 
and the formula (El), the following relations arise: 



Z = Sp{e- pH ) = J Db + (t)Db(t)e s , 



(^> = / Db+(t)Db(t) X A(0A(m)(fo+(f))' +1 (6+(t)) w ^(t))&r 1 We 5 /Z 

J 1,171=0 
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where the action S is defined by the form of the Hamiltonian H: 



r p ( f)h + (t) 1 
S = / dt(j2^^b t (t)--J2(X lk + Y tk )bt(t)S lk (t)b k (t) 

J ° \ i m 2 i+k 

- 7 E - Y lk )bt(t)bim 3 it) + h.c) 
^ i+k 

oo 

- £Z* £ a(l)a(m)(btmK(t)) m b\m(t) 

iy^k l,m=0 



1 ( 00 

- ~ E (h ix - ih iy ) £ + h.c. 

z ; ! 1=0 



- E^E«(0(^ + W)^W • 

i 1=0 / 

All other correlators can be obtained in the same manner and give us the bosonic 
functional integral representation which is free of constraints and limiting procedures 
The functional integral form then allows the simplest approach to the derivation of 
diagram technique rules which are standard ones for the problems in question. It 
is tempting to note that this technique is much less complicated and much more 
straightforward than the spin operator technique and is very natural for the con- 
sideration of problems concerning Bose-condensation (Long Range Order) in the 
system just using the standard Bogoliubov's approach to the subject. 
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